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ABSTRACT 

We report on a systeniatie and sensitive search for pulsars and transient sources in the nearby spiral 
galaxy M33, conducted at 1.4 GHz with the Arecibo telescope's seven-beam receiver system, ALFA. Data 
were searched for both periodic and aperiodic sources, up to 1000 pc cm"'^ in dispersion measure and on 
timescales from ~50 /xs to several seconds. The galaxy was sampled with 12 ALFA pointings, or 84 pixels 
in total, each of which was searched for 2-3 hr. We describe the observations, search methodologies and 
analysis strategies applicable to multibeam systems, and comment on the data quality and statistics of 
spurious events that arise due to radio frequency interference. While these searches have not led to any 
conclusive signals of periodic or transient nature that originate in the galaxy, they illustrate some of the 
underlying challenges and difficulties in such searches and the efficacy of simultaneous multiple beams 
in the analysis of search output. The implied limits are ;£ 5 /LtJy Mpc^ in luminosity (at 1400 MHz) for 
periodic sources in M33 with duty cycles s5%. For short-duration transient signals (with pulse widths 
S 100 /is), the limiting peak flux density is 100 mJy, which would correspond to a 5-ct detection of 
bright giant pulses (~20 kJy) from Crab-like pulsars if located at the distance of M33. We discuss the 
implications of our null results for various source populations within the galaxy and comment on the 
future prospects to conduct even more sensitive searches with the upcoming next-generation instruments 
including the Square Kilometer Array and its pathfinders. 

Subject headings: galaxies: individual (M33) - pulsars: general - pulsars: individual (Crab pulsar, PSR 
B0045+33, PSR J0628+0909) - surveys 



1. INTRODUCTION 

Numerous pulsar surveys over the past four decades 
have led to an imp ressive tally of ~2000 p ulsars known 
within our Galaxy (jManchester et all [20051) . The Magel- 
lanic clouds, located at a distanc e of ^50 kpc , are known 
to harbour 20 pulsars (Manchest er et ani2006[ ). which are 
the most distant pulsars currently known. Finding pul- 
sars in galaxies other than our own still remains a major 
challenge, as the achievable sensitivities of most currently 
operational radio telescopes are well below that required 
to detect even the most luminous pulsars in nearby galax- 
ies. Nonetheless, discoveries of such extragalactic pulsars 
will be scientifically rewarding on multiple counts: in ad- 
dition to providing valuable insights into the pulsar pop- 
ulations and birth rates within the host galaxies, such 
pulsars, through their dispersion measures (DMs) and ro- 
tation measures (RMs), have the potential to probe the 
baryon density and magnetic field in the intervening inter- 
galactic medium as well as the interstellar medium of the 
galax ies in which they reside (jCordes et al.ll200^ ICordesI 
I2007D . 

Of all the existing instruments, the Arecibo 305-m radio 
telescope is the most promising to conduct such searches 
given its superior sensitivity. However, its inherent limita- 
tion in sky coverage restricts the number of suitable tar- 
gets. Within the Local Group of galaxies, the Tri angulum 
galaxy M33, a spiral galax y located at 0.95 Mpc (|U et al.l 
l2009t iBonanos et all [20061 ) . is a highly promising target. 



Its location away from the Galactic plane (|5| ~ 30°) and 
the face-on orientation mean only moderate values of dis- 
persion measure are expected due to Galactic or extra- 
galactic material and therefore the search is unlikely to be 
limited by dispersion or scattering at frequencies above ~1 
GHz. Previous Arecibo searches bv lMcLaughlin fc CordesI 
(120031) were made at a frequency of 430 MHz and targeted 
giant-pulse like signals over a large DM range of 0-250 
pc cm" 



As highlighted by ICordes fc McLaughUnI (|2003[ ). stan- 
dard periodicity searches hold limited prospects for find- 
ing pulsars out to such large distances. However, alter- 
nate strategies, such as the detection of giant pulses from 
pros pective Crab-like pulsars in the galaxy, seem promis- 
ing (ICordes fc McLauehlinl l2003t iMcLaughUn fc CordesI 
120031 ICordes et al.l l200ir Giant pulses are sporadic 
bursts of emission with amplitudes reaching as large as 
two to four orders of magnitude brighter than those of 
regular pulses. For the Crab pulsar, they are known to 
reach amplitudes of ~160 kJy at 430 MHz (ICordes et al.l 
120041) and ~45 kJy at 1.4 GHz (iBhat et al.l [20081) . Such 
"giant-pulse" emission, once thought to be unique to 
the Crab, has now been detected from several more ob- 
jects, including millisecond pul s ars (MSPs) B1937-f-21 and 
J1823-3021A (ICognard et all 119961: [Knight et al. l I2005D 
and the LMC pulsar B0540-69 ([Johnston fc Romanil 
l2003l) . The estimated brightness temperatures are known 
to reach up to 10^^ K for sub-ns giant pulses from the 
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Crab (IHankins fc Eilekl I2007D an d lO^^K for the MSP 
B1937+21 (jSoglasnov et al.ll2n04D . 

One of the distinguishing characteristics of giant pulses 
is that their pulse amp litudes follow a power-law distri- 
bution (jLundgren et alJllQQSi) . However, no upper ampli- 
tude cutoffs have e ver been observed for any of the ob- 
jects s tudied so f ar ("Lundgre n et al.|[l995l : ISoglasnov et"aLl 
12004 iKnight et al. 2006), implying that exceedingly 
bright pulses are potentially detectable if observations 
are made over sufficiently long durations. Unlike period- 
icity searches, search strategies of this kind do not require 
unrealistically long integration times per pointing. How- 
ever, they rely heavily on the chance occurrence of such 
ultra-bright pulses during observations and consequently 
t heir expected rates of occurren ce . 

ICordes fc McLaughlinI (l2003l ) discussed the methodolo- 
gies and strategies relevant for searching for individual, 
dispersed pulses in fast-sampled spectrometer data. While 
originally developed for searching for giant pulse-like sig- 
nals in the data, these techniques are fairly generic and 
hence applicable for se arches of other targets that may 
emit transient signals. iMcLaughlin fc Corded (j2003fl re- 
port on systematic searches conducted with the Arecibo 
and Parkes telescopes toward several targets including 
M33 and the LMC. Subsequent applications of these tran- 
sient search techniques led to the discoveries of rotating 
radio transients (RRATs; McLa ughlin et a l. 2006) and a 
powerful burst of unknown origin (jLorimer et al.l I2007D 
and, most recently, s everal more RRAT- like objects in the 
PALFA surv ey data ("Pcneva et al.ll2Q09l ) and in the Parkes 
survey data (jKeane et al. 20^! 

While the searches for extragalactic pulsars by McLaugh- 
lin & Cordes (2003) did not result in any convincing detec- 
tions, they highlighted some of the underlying difficulties 
and challenges involved in such searches. In particular 
radio frequency interference (RFI) poses a major impedi- 
ment in such searches. Impulsive and powerful RFI bursts 
can potentially subdue, or even mask, the dispersed ra- 
dio pulses from pulsars, and ofte n tend to manifest as 
spurio us signals at non-zero DMs. IMcLaughlin fc CordesI 
(|2003[ ) found several tantalizing candidates in their data, 
however the single-pixel nature of their observations meant 
limited capability to distinguish them from spurious sig- 
nals of RFI origin, and therefore greatly limited their 
ability to conclusively establish a plausible a strophysical 
origin for any of them. In a follow-up work, I Bhat et al.l 
([200H) conducted exploratory searches toward M33 using 
simultaneous observations made with Arecibo and and the 
Green Bank Telescope (GBT). While the total observing 
time was rather too small and the sensitivity of the GBT 
too low to provide any detections, they demonstrated the 
efficacy of the dual-site aspect as a means of identifying 
and classifying various kinds of RFI (including that local 
to one or the other telescope and also satellite interference 
that is common to both). In any case, persistent RFI 
signals may adversely impact the detection sensitivities, 
as they tend to raise apparent detection thresholds, which 
m eans reduced s e nsiti vities to weaker signals. 

ILorimer et al.l (j2007i) reported a strong pulse near the 
direction of the Small Magellanic Cloud (SMC) which 
had a dispersion measure (DM) too high to originate 
from the Milky Wa y or fr om the SMC. More recently, 
iBurke-Spolaor et al.l (|201ll ) argued that it was swept- 



frequency RFI that convincingly approximated the ex- 
pected characteristics of an extragalactic radio transient. 
Regardless of such details and the nature of its origin, the 
report has stimulated great interest in probing the tran- 
sient radio sky on sub-second time scales. 

In § 2 we discuss the feasibility of using the Arecibo 
telescope to search for periodic and transient signals from 
M33. § 3 describes the data acquisition, § 4 the processing 
for transients and § 5 the periodicity search. We ana- 
lyze the results in terms of source populations in § 6 and 
summarize future possibilities in § 7. Our conclusions are 
presented in § 8. 

2. ARECIBO SEARCH SENSITIVITY AND SOURCE 
DETECTABILITY 

Arecibo's seven-pixel receiver system ALFA (Arecibo L- 
band Feed Array) offers an exciting opportunity to revisit 
such searches. Specifically, the combination of high sensi- 
tivity, large bandwidth and multiple beams provides both 
large instantaneous sky coverage as well as a powerful dis- 
crimination against spurious signals of RFI origin. Since 
2004, the system has been exploited for multiple large sur- 
veys i ncluding tho se for pulsars and sho rt-duration tran- 
sients CCordes et a l. 2006; Denc va~et al.| [2009). In this pa- 
per, we report on the first systematic searches conducted 
with the ALFA system toward the nearby spiral galaxy 
M33. Our searches targeted both periodic and aperiodic 
sources over a wide range of parameter space. Data quality 
is found to be remarkably good, with only a small fraction 
of the data corrupted by RFI. We discuss the efficacies of 
search strategies and methodologies that take advantage 
of ALFA's multiple beams and discuss the implications of 
our null results for various classes of source populations in 
the galaxy and place upper limits on their luminosities. 

Before we delve into the details of observations and 
data processing we briefly discuss different types of sources 
and signals that our searches are potentially sensitive to 
and comment on their detection prospects and the achiev- 
able sensitivities at a distance of ~1 Mpc. These in- 
clude: (i) ultra-luminous pulsars, (ii) bright giant pulses 
from Crab-like pulsars, (iii) rotating radio transients that 
emit ultra-bright bursts, (iv) transient radio emission 
from magnetar-like objects, detectable either as individual 
pulses or through periodic emission, and (v) other exotic 
bursts of Galactic or extragalactic origin. While our sur- 
vey was originally motivated by (i) and (ii), our processing 
techniques and strategies are fairly general and thus tar- 
get a much broader range of objects including the source 
classes of (iii) to (v). 

2.1. Periodicity Search 

For periodicity searches, the achievable sensitivity with 
ALFA is given by a modified form of the radiometer equa- 
tion: 
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where Lav is the period-averaged pseudo- luminosity at 1.4 
GHz (L band) in units of ^JyMpc^; given by Sav-D^, 
is the period-averaged flux density and D the 
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Fig. 1. — Sensitivity curves for periodicity searches with the 
ALFA system, plotted as minimum detectable pulsar luminosities 
at a distance of ~ 1 Mpc for a 8-cr detection threshold; the solid 
and dashed lines correspond to 1 hr and 2 hr integrations, respec- 
tively; with the thick (red) and thin (green) lines corresponding to 
the center and outer beams. Observing bandwidth assumed is 100 
MHz and the center frequency is 1.4 GHz. The filled circles de- 
note the 12 highest luminosity pulsars known within our Galaxy 
(i'l400~ 10'^ mjy kpc^ ), if they were located at the distance of M33 
(-Dm33 ~ 1 Mpc). 



distance; Dm33 is the distance to M33. We note that 
1 /iJy Mpc^ is equal to 10^ mJy kpc^, near the high end of 
pseudo-luminosities of Galactic pulsars at 1.4 GHz. Other 
quantities are: P, the rotation period; Wp , the mea- 
sured pulse width (thus Wp /P is the duty cycle); S/N, 
the signal-to-noise ratio; Sgys , the system sensitivity in 
Jy (given by T^ys / G where Tgys is the system tempera- 
ture and G is the telescope gain in K Jy~^ ); Av and Tint , 
the total bandwidth and integration time over which the 
data are recorded; and Npoi , the number of polarizations. 
The distance to M33 was previously estirnated to be ap- 
proximately 840 kpc (jFreedman et al.lll99lLl200lD. but has 
recently been revise d upward to 950 kpc ( U et al.ll2009t 
iBonanos et al.l[2006[ ). For th e nominal sensitivity parame- 
ters of the ALFA system fsee lCordes et al.l (|2006D l T^ys = 
30 K and G=10.4 KJy~^on-axis (i.e., center beam) and 
8.2 K Jy~^ off-axis (i.e., outer beams); hence, Sgys = 3.7 Jy 
and 3 Jy as the theoretical values for the outer and center 
beams respectively. 

The achievable sensitivities are shown in Figure [1] where 
the minimum detectable luminosity is plotted against 
the duty cycle of the pulsed emission, along with the 
~10 most luminous pulsars known within our Galaxy 
(Lav ~ 1 A* Jy Mpc^ ) . Of the 1880 pulsars known so far, 
nearly 75% have Lav measured at 1400 MHz. Interestingly, 
two pulsars - PSR J1644-4559 and PSR J0738-4042 with 
Lav in the 6-8 fiJy Mpc^ range - will be well above the pre- 
dicted detection curves (assuming 1 hr integration and 100 
MHz recording bandwidth) even when moved to a distance 
of 1 Mpc, with marginal detections possible for two others 
for longer (2 hr) integration. For an 8-a detection thresh- 
old, this corresponds to a minimum period-averaged flux 
density Sav ^ 4 /iJyin a 2 hr integration. So, if M33 har- 
bors a pulsar population with a similar luminosity distri- 
bution, indeed there should be at least a few pulsars that 
are potentially detectable with ALFA's sensitivity and in- 
tegration times of ~2 hr. 



2.2. Transient Search 

The search sensitivity for transient objects is given by a 
simpler form of the above radiometer equation, 

^ _ ( Dm33 \^ ( '5'pk \ / Ssj ^ 
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where S'pk is the peak flux density and Wp is the measured 
width of the detected pulse. This measured width is de- 
termined by the intrinsic pulse width, scatter broadening 
from multipath propagation, and dispersive smearing due 
to the instrument and any error in DM. Assuming Wint , 
Atscatt , Atdm.ch and Atdm, err denote these terms, the net 
pulse width measured is given by 
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Toward M33's direction (/ = 133.6°, h = -31.3°), the 
predicted Galactic co ntribution to th e dispersion measure 
is only ~50 pc cm~'^ (jCordes fc Lazio 2002) and therefore 
the pulse broadening is negligible at frequencies of 1-2 
GHz (Atscatt ~ 0.1 /is). Even if the model prediction is 
off by a factor of two, the dispersion smearing (Atdm,ch ) 
expected is ~100 /is for a frequency channel width of 
400 kHz. Thus, pulses that are intrinsically very nar- 
row (Wint ~ 100 /is) will be seen as '--^100 /is pulses, while 
broader ones will encounter little smearing from dispersion 
or scattering along the line of sight. 

Based on the above equation, a giant pulse that is as 
bright as 20 kJy at the Crab's distance (~2 kpc) will 
be potentially detectable at M33's distance (^1 Mpc) 
at the 5-cr level (i.e., a peak flux density, Spk ~100 
mJy). Thus, detection prospects critically depend on 
the rates of occurrence of such large-amplitude pulses. 
While several authors have recently investig ated the am- 
plitude distributions of Crab giant pulses ([Cprdes et all 
120041: IPODOV fc Stappers"2007': 'Bhat et al."2008'), the most 
extensive work is by Lundgr en et al . (1995) who analyzed 
70 hr of observation from the Green Bank 43-m telescope 
at 812 MHzQ Unfortunately, the distribution is poorly 
constrained at such large amplitudes due to small number 
statistics. In any case, within the statistical errors, the 
rate at which such bright pulses may occur is rather low 
(< 1 hr^^). Counter to this is the possibility of pulsars 
that emit giant pulses brighter than that of the Crab, as 
it would be surprising if the Crab turns out to be the one 
that emits the brightest giant pulses in the local Universe. 
Indeed, giant pulses with amplitudes as large as 50-100 
kJy have been detected at 1-2 GHz (or at least i n the 
L band) (|Popov fc Stappersl [20071: iBhat et ai]|2008[ ). but 
these are relatively narrow (;S 10 /isin duration) and so 
will be detected at a reduced sensitivity if the data are 
sampled at a rate of 50-100 /us, as is typically the case 
with pulsar search observations. 

Rotating radio transients are a newly discovered class 
of objects with amplitudes ranging from 10 mJ y to 4 Jy 
and pulse widths from 1 to 100 ms jM cLaughl hi et al.l 
l2006t iDeneva et ahl 120091: iBurke-Spolam fc Bailee! 120101) . 
While the ALFA receiver is almost an order of magnitude 



ICordes &: McLaughlinI II2003I ) presents a scaled version of their re- 
sults for a pulse width of 100 lis. 
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Fig. 2. — Sensitivity curves for short-duration transient searches 
with the ALFA system, plotted as minimum detectable peak flux 
density (Spk min) when the matched filter width is equal to the pulse 
width. The solid curve is for a 5-(t detection threshold and the sen- 
sitivity of the center beam; this curve will shift slightly upward or 
downward (i.e., dashed lines) for outer beams or for lower detection 
thresholds (e.g., 3-cr). Observing bandwidth assumed is 100 MHz 
and the center frequency is 1.4 GHz. The plotted symbols denote 
transient sources known within our Galaxy (i.e., giant-pulse emit- 
ters, radio magnetars and rotating radio transients), after scaling 
their measured peak flux densities to M33's distance. 



more sensitive than the Parkes multibeam (PMB) receiver 
(ICordes et aLll2006D . the PALFA survey uses much shorter 
dwell times than the PMB survey (67-268 s for the former 
vs. 2100 s for the latter). This difference might explain a 
general tendency where the PALFA survey detections are 
relatively weaker and tend to occur more frequently. The 
combination of ALFA's high sensitivity and our longer du- 
ration per pointing (2-3 hr) makes our search theoretically 
more sensitive to objects over a wider range of brightness 
and event rate. 

The achievable sensitivities for transient searches with 
ALFA (for a 5-a detection threshold) are shown in Fig. [2j 
along with the expected peak flux densities for currently 
known radio transients, if they are moved to a distance of 
Mpc. Evidently, the maximum achievable sensitivity 
is currently the fundamental limitation in finding such ob- 
jects at M33's distance. With even the brightest known 
objects well below the detection threshold, prospects are 
not good unless there is a population of much brighter 
objects in the galaxy. 

Another interesting class is magnetar-like objects that 
may turn "on" at radio wavelengths during their active 
phases. These are young, long-period pulsars with ultra- 
strong surface magnetic fields (~ 10^'* G) whose emission 
is thought to be powered by the r earrangement and de- 
cay o f their extreme magnetic fields (jPuncan fc Thompsonl 
Il992l ). They provide crucial links between normal pul- 
sars and the clas s of anomalous X-ray pulsars (AXPs; 
iMereghett] l|2008M Their radio fluxes ten d to vary dra- 
matic ally over timescales of days or weeks (jCamilo et al.l 
|2006() and thus they can potentially be detected through 
either periodicity searches or through searches of single 
pulses, depending on the brightness. With just three ob- 
jec ts known so f a r (XTE J18 10-197 and IE 1547-5408 
bv ICamilo et al.l (l2006l 1200711 and PSR J1622-4950 by 
iLevin et al. ( 2010f )). their population remains poorly con- 
strained. However, if we take these three objects as repre- 
sentative of the class, their detection prospects are rather 
poor at M33's distance. The brightest of the three ob- 



jects has Lav, 1400 - 0.15 /iJyMpc^ , i.e., almost two orders 
of magnitude below the detection limit at a distance of 
'^l Mpc (see Fig. [T]). The brightest pulse detected from 
XTE J1810-197 has S'pk,i4oo^ 40 Jy, Wp S 10 ms, which 
is still an order of magnitude smaller than that required 
to enable a detection at '^1 Mpc (Fig. [5]). 

In addition to the above-described source classes, there 
are likely to be several as-yet-unknown types of sources 
that emit bursts of radio emission over short time dura- 
tions. Among such hypothetical source classes are anni- 
hilating black holes and gravit ational wave events (jReesI 
Il977t [Hansen fc Lvutikovii200l[ ). The recent detection of 
a power ful burst (with a flu x of 30 Jy and a width of 5 
ms) by iLorimer et al.l (j2007f) triggered much interest in 
such exotic sources. The non-detection of more such sig- 
nals in the PALFA survey argued again st an astrophysical 
nature of its or igin (Dc neva et al.|[2009[l and more recently 
iBurke-Spolaor c t al. ( 20111 ) convincinglv described it to be 
of terrestrial origin. In any case, if any such bursts were 
to occur during our observations, they will be easily de- 
tectable. 

3. SURVEY POINTINGS AND OBSERVATION PARAMETERS 

Observations were made during September-December 
2006 using Arecibo's ALFA system which allows simul- 
taneous data collection from seven fields, each 3.5' 
(FWHM) across. Taking into account the hexagonal ar- 
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Fig. 3. — ALFA pointings used for the M33 search overlaid on an 
optical image of the galaxy from the Sloan Digital Sky Survey (York 
et al. 2000). The optical radius of the galaxy is roughly 30 arcmin. 
There are 137 supernova remnant candidates in the galaxy which are 
uniformly distributed out to a distance of 15 arcmin from the center 
(Long et al. 2010). As discussed in the text, the optimal sampling 
possible with a maximum of 12 ALFA pointings still leaves a small 
part of the galaxy sampled with a gain less than half-power of the 
on-axis gain. The three closely spaced ALFA pointings, which fill 
in the area covered by the ALFAs beam pattern as projected onto 
the sky (cf. Cordes et al. 2006 for details), are denoted as A, B, 
and C in Table 1. The galaxy coverage from our 12 ALFA pointings 
almost fully overlaps with that of the 16 pointings of McLaughlin 
&; Cordes (2003) used for their 430 MHz single-pixel search. 



Arecibo Search 



5 



rangement of the beams on the sky and the near sidelobes, 
the combine d power pattern is approximately 24' x 26' 
(jCordes et al. 2006). Thus, in order to cover a given field 
of the sky (with a sensitivity better than the half-power 
point of on-axis gain), three closely spaced ALFA point- 
ings are needed in order to fill in the holes in the beam 
pattern projected onto the sky (see Fig. 3 of ICordes et al.l 
(|2Q06)). 

Given M33's extent of roughly 73' x 45', 12 such ALFA 
pointings are required to achieve a dense sampling of the 
galaxy. Details of pointings and observation parameters 
are summarized in Table [TJ Fig. [3] shows the tiling scheme 
and pointings used for our survey. While this represents 
an optimal sampling possible with 12 ALFA pointings, it 
still leaves a small fraction of the galaxy unsampled. A 
complete sampling as per the ideal tiling scheme as used 
in the PALFA survey would require many more pointings, 
with a significant number of beams pointed away from the 
galaxy, and hence is not an efficient use of telescope time. 
We therefore optimize our search strategy by limiting the 
search to 12 pointings, which cover approximately 95% of 
the galaxy. 

For eight of these pointings, data were collected for 
2 x 1 hr durations, while four other pointings were ob- 
served for 3 X 1 hr. Data were recorded using the Wide- 
band Arecibo Pu lsar Processor (WAPP) spectrometers 
(jDowd et al.ll2000D . which were used to synthesize 256-lag 
autocorrelation functions for both three-level quantized 
polarization channels; correlation functions for the two 
channels were summed before writing out to disk. We ob- 
served a 100 MHz passband centered at 1440 MHz in each 
of the seven telescope beams; therefore, each frequency 
channel is approximately 390 kHz in width. In order to 
accommodate the constraints of the available data storage 
resources, approximately half of these data were sampled 
at a rate of once every 64 /is while the remaining half at a 
slightly slower sampling of once every 100 us (see Table[T]). 
This resulted in an aggregate data rate of approximately 7 
X 30 GB per every hour of observation, amounting to a to- 
tal data volume of 6 TB from 15 observing sessions. These 
data were transported to the Swinburne supercomputing 
facility where the search processing was carried out. 

The RFI environment at the Arecibo Observatory is 
known to be a strong function of time of the day, with a 
tendency for minimal occupancy in time and frequency be- 
tween midnight to dawn hours. Fortunately, all our obser- 
vations were made during such relatively radio-quiet hours 
of the day when RFI from terrestrial sources is generally 
at a minimum. As we highlight in a later section, overall 
data quality is found to be excellent, with only a small 
fraction (sO.1%) of the data affected by RFI. Moreover, 
the scheduling of our observing sessions were coordinated 
with those of the PALFA survey, usually centered at 19 hr 
LST (i.e., a few hours prior to our sessions). Successful 
detections of test pulsars at the beginning of our observ- 
ing sessions and the discovery /confirmation observations 
of several new pulsars made from the PALFA survey data 
collected during this period thus readily provide important 
integrity checks of the full observing system and processing 
pipelines. 



4. SEARCHES FOR TRANSIENT SIGNALS 

As outlined in §2, our search processing targeted two 
kinds of signals: (i) bursts of emission from giant-pulse 
emitters such as the Grab, or intermittent sources such 
as RRATs or magnetars, or one-off bursts from hitherto 
unknown sources (e.g., Lorimer et al.); (ii) periodic sig- 
nals from highly luminous pulsars or active magnetars. 
For searching for transient signals of the first category, 
the procedu res that we adopted are very similar to those 
described in iDeneva et al.l (|2009t) . except that our search 
extends to a wider range in parameter space as the time 
scales we search for range from 64 /xs to several seconds. In 
summary, raw data were dedispersed over a wide range in 
DM and each dedispersed time series is searched for indi- 
vidual pulses with two different time domain algorithms: 
(i) matched filtering, which is the most commonly used 
technique in single pulse searching, and (ii) time domain 
clustering, which employs an algorithm similar to friends- 
of-friends and thus sensitive to clustering of samples above 
threshold in time. 

In the matched filtering technique, a pulse template is ef- 
fectively convolved with the dedispersed time series, where 
a template is chosen to reflect the diversity of pulse shapes 
that are searched. For example, unscattered pulse shapes 
can be modeled as either a single or a sum of two or more 
Gaussians whereas asymmetric pulse templates would be 
a better approximation for scattered pulses which tend 
to have exponential tails. In practice, we approximate 
matched filtering by progressive smoothing of the time se- 
ries by adding up to 2" neighboring samples and select- 
ing events after a set threshold after each iteration. In 
light of the marginal detection prospects anticipated in our 
searches, we use a fairly low threshold of 3.5 tr for event 
detection. As the smoothing is done in pairs of samples, 
the pulse template used is effectively a boxcar of length 2" 
samples. 

The friends-of-friends method relies on a completely dif- 
ferent approach, where the dedispersed time series is pro- 
cessed sequentially and if an event above the threshold is 
found it is designated as the first of a cluster. A cluster 
of events is then augmented while successive samples are 
found to be above the threshold. The brightest sample of 
a cluster is recorded as the event amplitude and the total 
number of samples in the cluster as its width. Such an 
approach has a maximum effective time it is sensitive to, 
which is related to the block size of the data used in calcu- 
lating statistics of the time series. In our case, it is limited 
to 4096 time samples while searching for short-duration 
events (;£ 1 s); i.e., 262-410 ms depending on whether the 
sample time used in the observation is 64 or 100 /its . 

4.1. Data processing 

We dedisperse the raw data over a wide range of 0-1000 
pc cm~'^ in DM using 1272 trial values that were optimally 
chosen to span this range. Although the DM to M33 is es- 
timated to be only ^50 pc cm~^, we set an upper limit 
of 1000 pc cm~^ in order to account for any pulsars that 
might be in or behind H II regions within M33 and to 
account for any potential extragalactic pulses originating 
from behind M33. For searches based on the matched fil- 
tering scheme, we used boxcars of widths 2", where n rang- 
ing from to 17, thus covering a large range in timescale 
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Fig. 4. — Diagnostic plots from transient search pr ocessing: (Left) Detections of individual pulses from J0628+0909, a RRAT-likc object 
discovered in the PALFA survey llCordes et al.l [20061) : (Right) Detections of single pulses from a known pulsar PSR B0045+33 that was 
routinely observed as a test pulsar at the start of every observing session. In each set, the top panels (left to right) show the statistics of 
events above 5-ij against S/N and DM, the scatter plot of S/Ns and DMs, and the mapping from the DM channel to DM value. The bottom 
panel shows the summary of events (S/N > 4) in the time-DM plane. Multiple bright pulses are seen at the s periodicity of the pulsar 

B0045+33; the brighter pulses tend to peak near the pulsar's true DM fa 40 pc cm~'^ (DM channel = 15). As the pulse width is 3 X larger 
than the dispersion sweep across the 100 MHz band, most pulses tend to appear at a large range in DM. The signal strength tends to peak 



near the true DM of the object, 



pc cm ^ , for J0628+0909. The size of the open circle is proportional to S/N. 



from 64 us to 8.4 s. Even though the detected emission 
from fast transients known so far have timescales S 100 
ms, this extended range allows a much more comprehen- 
sive search, as it targets even unusually broad pulses or 
strongly scattered pulse shapes. This strategy based on 
the use of symmetric templates is however non-optimal 
for scattered pulses (which are asymmetric due to their 
exponential tails), which will therefore be detected with a 
reduced sensitivity. 

For processing efhciency, the raw data were broken into 
short blocks of 135 seconds each, which is set by the 2 
GB file size limitation of the data recording. Moreover, 
the matched filtering procedure was done in two passes: 
the first pass with n ranging from to 7 (i.e., pulse widths 
ranging from 64/100 /iS to 8.192/12.8 ms depending on the 
time sampling), and the second from n—8 to 17 (i.e., pulse 
widths ranging from ^10 ms to ^10 s), thus spanning 
almost five orders of magnitude in timescales. The com- 
bination of short data blocks and two separate passes for 
short and long events (s 10 ms and ^10 ms, respectively) 
avoids overcrowded diagnostic plots and thus a much more 
efficient process of sifting through the s earch outputs. 

As highlighted bv lDeneva et al.l (j2009r ). the time-domain 
clustering method is complementary to the matched fil- 
tering technique, with the major advantage being its abil- 
ity to efficiently suppress many false positives that arise 
from RFI sources. This is because a broad RFI burst will 
most likely be detected as a single event by this method 
as opposed to a large number of events as is the case with 
matched filter method. A limiting factor however is that 
the sensitivity is limited by the sampling resolution as 
there is no progressive smoothing of the time series and 
therefore is it less sensitive to weak, narrow pulses. Any 
genuine signals of astrophysical origin in the data will most 
likely be detected by both the methods, albeit at different 
sensitivities depending on the pulse width. This is the case 
in the PALFA survey. 



4.2. Analysis of Search Outputs 

Example diagnostic plots from our search processing are 
shown in Fig. [H These show all candidate events with 
signal-to-noise ratio (S/N) larger than the 5-cr threshold 
and their statistics in the form of distributions of S/N 
and DM as well as a scatter plot of S/N vs. DM. Any 
genuine signals of astronomical nature will have distinct 
signatures on such plots; for instance, they will be de- 
tected over a contiguous range of DMs, peaking near the 
true DM and with a signal strength that falls off smoothly 
with an increasing departure from the true DM. Ideally, in 
the absence of RFI, there will be a clear peak in the DM 
distribution as well as in the scatter plot (S/N against 
DM), peaking near the true DM of the signal. The width 
of the peak vs. DM is broader for broader pulses. The 
density of events on the scatter plot will depend on the 
number of similar events detected in the data. In con- 
trast, RFI with narrow time structure will peak at zero 
DM and fall off rapidly at larger DMs unless it has intrin- 
sic swept-frequency structure. If very strong or temporally 
wide, the RFI can appear over the full range of DMs. Such 
strong and impulsive RFI bursts can potentially modify or 
even mask the signatures of any real signals present in the 
data. 

In principle, a careful assessment of these diagnostic 
plots should allow identification of genuine signals vs. those 
arising from RFI; however, the efficiency of such a strategy 
critically depends on the nature and strength of RFI. As is 
well known, impulsive RFI can often mimic the signatures 
of real signals and may manifest as islands in the time- 
DM plot or as distributions which peak at certain DMs. 
Indeed, as we discuss in § 14.31 a more rigorous scrutiny 
that involves inspection of relevant raw data segments in 
the time-frequency plane in order to confirm the expected 
dispersion sweep will help detect and eliminate all such 
cases with high enough signal-to-noise for frequency-time 
analysis. However, the large number of events typically 



Arecibo Search 



7 



s/n > ft: 

L't'* 9 ..■'".Ax 




C3 



E 




Li) iDoo 1^.* 



500 H'^a? 





eO IIHD 

E 

P 500 

pq 



1 




■ ■■- ""ip 





10 IP' ifl' 10 *o ao 40 
1, n/'^ 



J_' " 1 " " 1 " 


'1""!; 


1 ■ >■ ■!■ ■■>): 


J 1 1 ■ 1 1 ' 1 1 1 1 1 


.|. .M|. 






1' " r " ' 1 " ■ ' 1 " " t 










iH.|mi|ii 








1, ., 1: 




■i 











1 e 3 



Fig. 5. — Example diagnostic plots from event analysis and classification: An isolated event from observations made on MJD 53992 toward 
the pointing M33.4.C, detected with a signal strength (S/N ~30) that peaks near DM~200 pc cm""^ in beam 5, indicating a broad, dispersed 
pulse. A closer scrutiny reveals this tantalizing event to be of instrumental origin, largely a manifestation of an abrupt increase in power level 
that lasted for a short duration. The localization in DM is intriguing and an explanation is given in § 4.4. Each row shows results from one 
ALFA beam along with dedispersed time series. From left to right, panels show: events with S/N > 5 vs. DM channel and time, number 
of events vs. DM channel, event S/N vs. DM channel, mapping from the DM channel to DM value, and time series of a short data segment 
around the event. 



detected in our searches make such an approach both te- 
dious and impractical. 

Fortunately, multiple beams of the ALFA system can 
be well exploited for this purpose given the significant im- 
provement it offers for event analysis. However, such a 
filtering scheme based on multiple beams is not straight- 
forward because scattering off telescope structures can in- 
troduce RFI into the telescope optics. Depending on the 
telescope orientation, RFI bursts may be detected in all, 
some or none of the ALFA beams. We therefore adopt a 
scheme whereby search outputs from multiple beams are 
compared and inspected visually, as described in the next 
section. 

4.3. Event classification and Scrutiny 

In order to aid a meaningful classification of events and 
their final scrutiny, we assemble various search output 
plots as described in § 14.21 for all 7 beams to generate 
multipanel plots as shown in Figs. [5] and |6l Such plots 
readily provide a comprehensive summary of all detected 



events in a given data set as well as data quality, and thus 
drastically improve the efficiency of event analysis. Several 
thousands of plots generated in this manner are included in 
a MySQL database which was later accessed and queried 
using a python-based interactive viewer for ranking and 
classification. We adopted a scheme whereby each data 
segment is first assessed and ranked, with ranks of 1 and 

2 (Class I and II) to reflect prospective candidates that 
require further careful scrutiny, while data that are devoid 
of any significant events (i.e., no clustering in the time- 
DM plane or in the statistics plots) were given a rank of 

3 and those which identified to contain at least one event 
of RFI origin were given a rank 4. The best candidates 
are thus of Class I, while those of Class II are considered 
promising. Results of this assessment are summarized in 
Table [21 The RFI bursts are found to typically last over 
time durations of a few seconds or less, and no more than 
a few instances detected in a given data segment. Thus 
only a small fraction of the data are infiuenced by RFI and 
therefore overall data quality is pleasingly good. 
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Fig. 6. — Example diagnostic plots from event analysis and classification: a spurious event from observations made on MJD 53995 toward 
the pointing M33.2.B, which appears at a range of DMs and signal strengths simultaneously through the center beam and the beams 2 and 
5. The event was classified to be of RFI origin, given a highly asymmetric nature of the S/N vs. DM curve (in beam 0), and its appearance 
(with much lower S/Ns) in two well-separated outer beams (2 and 5). Each row shows results from one ALFA beam along with dedispersed 
time series. From left to right, panels show: events with S/N > 5 vs. DM channel and time, number of events vs. DM channel, event S/N vs. 
DM channel, mapping from the DM channel to DM value, and time series of a short data segment around the event. 



Candidate events of Class I and II are then subjected to 
a more rigorous scrutiny, where short raw data segments (4 
seconds) centered at the time of the event were extracted 
for ah seven beams and displayed in the frequency-time 
plane (i.e., dynamic spectrum) in order to check the ex- 
pected dispersion sweep across the frequency band. This 
provides a crucial diagnostic for celestial vs. terrestrial ori- 
gin. These dynamic spectra plots are examined along with 
the dedispersed time series for the event's DM, to facili- 
tate a more comprehensive scrutiny of the candidate event. 
For the beam corresponding to the candidate event, fur- 
ther scrutiny is carried out by also inspecting the dedis- 
persed time series at two nearby DMs and at zero DM. 
This analysis is typically done at the original resolution of 
the data (64 or 100 /is). However, for candidate events 
with S/N:S 10, and thus too weak to be discernible in the 
frequency-time plane, data are suitably smoothed in time 
and/or frequency and then re-inspected. Such a procedure 
allows a rigorous scrutiny and the final assessment is done 
in this manner for all 169 candidate events that emerged 
from our search. 



Example plots illustrating this scheme are shown in 
Fig. [71 which shows the signatures of a real signal (an 
individual pulse detected from J0628-f 0909, a RRAT-like 
ob.iect that was orig inally discovered in the PALFA survey 
(jCordes et al.ll2006ti ') and an RFI event. 

In summary, our scheme for event analysis, which relies 
on multiple diagnostic plots, detailed comparisons with 
results from other beams and careful examinations of the 
raw data segments, have allowed us to conduct a thor- 
ough and comprehensive assessment of all prospective can- 
didate events. Unfortunately, we have not found any signal 
for which an astrophysical origin can be convincingly at- 
tributed. Instead, all such events (Class I and II) turned 
out to be intriguing manifestations of different kinds of 
RFI, with the pulse widths ranging from ~100 /isto sec- 
onds and DMs from ^10 to ^300 pc cm"'^ . The pulse 
shapes vary from highly spiky to rectangular, and from 
smoothly symmetric to highly exponential, or sometimes 
seen as trains of short bursts. More details on the na- 
ture of RFI and some quantitative analysis are described 
in §1121 
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4.4. Radio Frequency Interference and Data Quality 

Spurious events that originate from impulsive RFI are a 
major problem in transient detection. However, as high- 
lighted in previous sections, the overall data quality is 
found to be pleasingly good in our observations and only a 
small fraction of the data (sO.1%) are corrupted by RFI. 
Moreover, our extensive scheme of event scrutiny that re- 
lies on multiple diagnostics including the event signatures 
in the time-DM plane, signal strength vs. DM and its ap- 
pearance across multiple beams, gives us a powerful han- 
dle in the identification and elimination of the bulk of such 
RFI-generated events. In particular, the infamous, and of- 
ten periodic, strong RFI pulses from the Federal Aviation 
Administration (FAA) radar at San Juan airport, which 
have a 12-sec o nd pe riodicity and are ~1 second wide (see 
iDeneva et al.l (|2009[ ) for details), are practically absent in 
our data. This is most likely to do with the fact that our 
observations were carried out during the radio-quiet hours 
between midnight and dawn, when the FAA radar tends to 
be relatively less active. Furthermore, in the small number 
of instances when such radar pulses were seen, their nature 
is found to be rather aperiodic. Given such low presence 
of radar RFI, a rigorous excision scheme as discussed by 
iDeneva et al.l ([2009) is not warranted for our data. 

Table 2 provides a short summary of the number of can- 
didate events identified in our analysis and those which can 
be readily identified as RFI. Our assessment of an event to 
be of RFI origin is based on one or more criteria, such as 
a high signal strength at zero DM and its appearance in 
many or all DMs, as well as in several or all 7 beams. As 
seen from this summary, the total number of RFI events 
detected depends on the detection method, with the time 
domain clustering method comparatively more efficient in 
their detections. On the other hand, the matched filtering 
technique results in a much larger number (four times) of 
"interesting" events. In general, with any broad bursts 
(RFI or real), we would expect the matched filtering tech- 
nique to result in a larger number of events than the time 
domain clustering method (cf. § 14. ip , however the number 
of events will depend on factors such as the pulse width, 
block length and the number of boxcars used in process- 
ing. As discussed in § 14.31 all such promising candidates 
are subjected to a thorough scrutiny, involving the extrac- 
tion of the relevant raw data segments and their inspection 
both in the time-frequency plane as well as in time series 
at multiple DMs, across all seven beams (Fig. [7]). 

Figures [5] and [5] show representative examples of such 
interesting candidate events seen in our analysis. In many 
cases, the most striking feature is the detection over a 
contiguous range of DMs, with the signal strength peak- 
ing at a certain DM, and its occurrence in only one or a 
few beams. While at first glance such signatures may sug- 
gest a likely non-terrestrial origin, our close scrutiny re- 
veals that they are often manifestations of intriguing RFI 
pulses. For instance, broad, asymmetric pulses such as 
those in Fig. [3 with a significant rise time (~100 ms) and 
a pulse width (^250 ms) that is comparable to the block 
length used in our analysis, can potentially mimic a local- 
ization in the DM space, even though it is intrinsically at 
DM=0. The same may apply to jumps in power levels that 
last over durations of the order of block length or longer 
(e.g., Fig.[5])- The localization in DM may occur when the 



pulse width (or the duration of the event) is comparable 
to the block length used in processing and if at least one 
of the pulse edges occurs near the edge of the block. In 
such cases, de-dispersion results in smearing of the burst 
and, for certain DM values, the burst may be split into 
two successive blocks after de-dispersion. At both low and 
high DMs, much of the burst will still lie within a single 
block, which means larger apparent detection thresholds. 
Consequently, the burst may not be detected as an event. 
The localization in the DM space that arise in this man- 
ner is thus purely a processing artifact. This is especially 
the case in our processing pass targeted at "short" events 
(time durations S 10 ms) where the block length used in 
determining the noise statistics is restricted to 4096 time 
samples (~400 ms), which is comparable to maximal dis- 
persion sweep (278 ms at DM=1000 pc cm~'^ ). The extent 
of DM localization will primarily depend on the rise time of 
the pulse and, to a certain extent, other signal characteris- 
tics such as the pulse amplitude and width, and in general 
longer the rise time broader the extent of localization. We 
have inspected all such candidate events detected in our 
analysis and found none that may warrant attributing an 
astrophysical origin. In many cases they were recognized 
as RFI pulses from the FAA radar (sporadic or aperiodic 
in nature), and in some cases sudden jumps in intensity 
levels that are likely of instrumental origin. 

5. SEARCHES FOR PERIODIC SIGNALS 

Data were also searched for periodic signals, for which 
we used the PRESTO search codc0, following procedures 
similar to those adopted for the offline processing of the 
PALFA survey data. In summary, the raw data were first 
dedispersed over a large range of DM, 0-1000 pc cm~^ , 
subdivided into four groups, with successive ranges us- 
ing a larger step size and coarser time resolution. Specif- 
ically, a total of 1272 trial DMs were used for the search: 
with 624, 240, 264 and 144 trial DMs spanning the 0- 
1000 pccm"'^ range in steps of 0.3, 0.5, 1.0 and 3.0 
pc cm~'^ respectively. For the last three DM groups, the 
data were progressively smoothed and decimated (in time) 
by a factor 2, 4 and 8 respectively. 

A Fast Fourier Transform was performed on the result- 
ing dedispersed time series and the power spectrum was 
searched for periodic signals by summing up to 16 harmon- 
ics and applying a harmonic sum threshold of 6a. With 
our 1 hr long integration (> 2^^ samples) for each point- 
ing, this translates to processing 7 x 60 GB data per point- 
ing, typically producing hundreds of candidates per each 
beam. The large number of candidates that emerged from 
this search (11,879 in total) were loaded into a MySQL 
database and subjected to a careful inspection for any 
prospective candidate using an interactive viewer tool. 

A thorough scrutiny of these candidates revealed ^10% 
to be above an 8-a detection threshold as discussed in ij |2.1l 
Not surprisingly, a large fraction of them turned out to be 
of spurious origin, essentially manifestations of strong RFI 
pulses originating from the FAA radar (see § 14. 4p . Such 
candidates are generally detected with a roughly 6-s pe- 
riodicity or its harmonics (i.e., P ~ 3, 2 or 1 s) at low 
or moderate DMs (;S 50 pc cm""^), and inspection of the 
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Fig. 7. — Example diagnostic plots from event analysis: the top two rows show a spurious event at DM~70 pc cm~^ produced by impulsive 
RFI de tected through the center beam of ALFA, while the bottom two rows show a real signal from J0628+0909 discovered in the PALFA 
survey llCordes et al.|[2006l'l . In each set, the top rows of multiple panels show the dedispersed time series and dynamic spectra plots for all 
7 beams, and the bottom panel additional diagnostic plots for the beam in which the candidate event was detected. See the text (§ 14.311 
for further description. The time-frequency plot of the J0628+0909 pulse clearly displays frequency structure that is caused by interstellar 
scintillation. 
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Fig. 8. — Plot of pseudo-luminosity Si4ooD'^ against the spin pe- 
riod, assuming M33 harbors a pulsar population that is very similar 
to that o f the Milky Way. The da ta points are from the ATNF pulsar 
catalog dManchester et ani2005| ). scaled for a distance of 0.95 Mpc, 
where M33 is located. The solid curve is the 6(t detection threshold 
at 1.4 GHz based on the observing parameters described in the text 
and for identification of a preliminary set of periodicity candidates. 
The other twro curves indicate the detection threshold for a survey 
with the Square Kilometer array for a 1-hr integration at 2 GHz 
with a 1 GHz bandwidth (dashed line) and at 0.43 GHz with a 
0.1 GHz bandwidth (dotted line). A sensitivity A/T = 10* m^ K^^ 
is assumed for the SKA. 



relevant time series show muhiple instances of such RFI 
bursts and a folded pulse profile that has a characteris- 
tic asymmetric (exponential) shape. Only a tiny fraction 
(^2%) of the low S/N candidates (6-8 a) was found to be 
caused by the radar RFI. 

Of the remaining candidates, roughly 40% were found 
to be within a narrow region of the parameter space (DM 
S 5 pc cm~^ , P ;S 5 ms), while the majority (~60%) were 
weak, short-period candidates spanning a large range in 
DM. We inspected the diagnostic plots for all candidates 
at DM > and above 6-cr, and found none that satisfied 
the standard assessment criteria to be considered real (e.g., 
persistence of signal vs. time and vs. frequency, a clear 
peak in the signal strength vs. DM and in the period- 
period derivative search space, etc.). As seen from Fig.[Tl 
our integration times of 1-2 hr will translate to limiting 
flux densities of ^3-10 /iJy at M33's distance, depending 
on the duty cycle (Wp /P -1-10%). 

6. IMPLICATIONS FOR SOURCE LUMINOSITIES AND 
POPULATIONS 

Figure |S] shows the sensitivity curve for our analysis. 
The solid curve takes into account radiometer noise, and 
propagation effects in the Milky Way (disp er sion and scat - 
tering) using the NE2001 model (jCordes fc Lazid I2002D . 
At the Galactic latitude of M33 {b « —31°) and for our 
observation frequency (1.4 GHz), pulse broadening from 
residual dispersion smearing and from scattering in the 



Milky Way or from M33 are negligible. For simplicity, we 
have assumed that M33 harbors a pulsar population sim- 
ilar to that in the Milky Way. Only a few known pulsars 
are just above the detection curve, implying that if there 
were —2000 or more active pulsars that beam toward us 
from M33, we should have detected a few of them. How- 
ever, the mass of M33 is about 10% of the Milky Way's 
(e.g., Sakamoto ct al. 2003), so there may be fewer than 
2000 detectable pulsars in M33. A counter argument may 
be in order if the expected neutron star population in M33 
is larger. For instance, while we may expect the star for- 
mation rate in M33 to be roughly 10% of that in the Milky 
Way bas ed on mass argum ents, the measured rate is about 
10-25% (|Abdo et al.ll2010[ ). suggesting that there may be 
more neutron stars than we may expect based on the 10% 
rate. 

The figure also gives curves for the Square Kilometer 
Array at frequencies of 0.43 and 2 GHz. We have assumed 
a sensitivity of 0.28 Jy, corresponding to an effective area 
divided by system temperature of lO'* m^ K~^. The nomi- 
nal sensitivity for SKA Phase 1 is appr oximately the same 
as that of Arecibo ([Garrett et al.ll2010t) . while the full SKA 
may have a sensitivity that is several times to a factor of 
ten larger than that of Arecibo. 

A simple upper bound can be placed on transient 
sources in M33. We estimate the volume sampled per 
telescope beam diameter 0b ~ 3.5 arc min as AVi = 
(7r/4)(Z?M33^b)'i?eff « 0.08(i?eff/0.1 kpc) kpc3, where 
Hes is the extent of the galaxy sampled along the beam 
direction (assuming a face-on orientation of the galaxy). 
With a total aggregate search time of 30 hr (c.f. Ta- 
ble 1), the limit on the event rate per unit volume is 

(7AFiTtotai)"^ « 0.06(0.1 kpc/i?eff) events kpc'^ hr'^ 
for events that satisfy S'pk~23.5 mJy Wms^^- The Crab 
pulsar emits pulses with amplitudes consistent with this 
expression (when instrumentally broadened to 100 fis and 
accounting for the ratio of distances of the Crab pulsar 
and M33) about once per 10 hr, so the event rate limit cor- 
responds to an upper bound of about one Crab-pulsar-like 
object per kpc^ in M33. 

It is well known that pulsars have high space veloci- 
ties with a distribution that extends w ell beyond the local 
Gala ctic escape velocity — 500 km s~^ (jArzoumanian et alj 
120021) . For the 84 beam areas of our survey, the volume 
out to a distance of 5 kpc (i.e., the distance traversed by 
a canonical pulsar moving at 500 km s~^ in its 10 Myr 
lifetime) above the Galactic disk is ~ 3 kpc^. Assuming 
roughly 50% of the known — 2 x lO'^ active radio pul- 
sars that beam toward us are on escape trajectories, there 
should be about three pulsars in the search volume. The 
lack of any detections is consistent with this estimate given 
that most of the pulsars will not reach 5 kpc during their 
radio-emitting lifetimes. 

Lorimer et al (2007) reported a single 30- Jy, 5-ms event 
from their survey of the Magellanic Clouds that covered 
—9 dcg^ at high Galactic latitudes. Assuming a source 
population that is homogeneously distributed in Euclidean 
space, Deneva et al. (2009) estimate —600 such events 
above the detection threshold of the PALFA survey. These 
limits follow from the calculations where the number of 
detected events N (x fl Tobs Dmax (where is the solid 
angle coverage of the beam and Tobs is the survey dura- 



12 



Bhat et al. 



tion) and by ignoring the propagation effects that may 
hmit the detectability of bursts from distant sources. The 
maximum distance to which an event of a given peak flux 
density S'pk and the pulse width W is detectable is given by 

Dmax oc (m Ssys )'^^^Aj/ "■^/^ (cf. Equation [2l where m is 
the detection threshold in units of a). The sensitivities of 
our single-pulse searches are similar to those of the PALFA 
survey, however our total survey duration is substantially 
smaller (30 hr vs. 461 hr). We thus estimate '^40 pulses 
above our detection threshold of '^20 mJy (cf. ? 12. 2p . In 
any case, the non-detection of such plausible extragalactic 
events in our survey reinforces the conclusion of Deneva 
et al., effectively precluding a source population that has 
a homogenous and isotropic distribution in the sky. 

7. FUTURE PROSPECTS 

Detection of extragalactic pulsars and transients is in- 
evitable with continued observations with existing tele- 
scopes and observations with forthcoming instruments. 
As discussed in [JTl the Crab pulsar's giant pulse ampli- 
tude distribution is a power law, with no upper limit hav- 
ing yet been detected. Further, the rate of giant pulses 
from the Crab pulsar is also a power law in pulse en- 
ergy, with brighter pulses being less frequent than fainter 
pulses. However, because no upper limit to the pulse am- 
plitude distribution has been observed, longer dwell times 
imply a higher probability of brighter, and therefore more 
detectable, pulses. Further, it is plausible that objects 
brighter than the Crab pulsar exist. 

There are now a suite of centimeter- and meter- 
wavelength facilities either in commissioning or under con- 
struction that suggest optimism for the future detection 
of extragalactic pulsars and radio transients. These in- 
struments will provide higher raw sensitivities (i.e., lower 
Ssys) and/or larger bandwidths than we have had available 
in this Arecibo search. For specificity here, we focus on 
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D (Mpc) 

Fig. 9. — Distance to which a bright giant pulse — having a flux 
density of 20 kjy, from the Crab pulsar (2 kpc) — could be detected 
with future instruments. See also Equation (2) and Table 3. The 
vertical dashed line indicates the distance of M33. For clarity, the 
curve for Arecibo is not labeled, but it is the purple curve between 
the curves for FAST and SK A- 1 . The dashed purple curve is for our 
current Arecibo search. We do not plot above a signal-to-noise ratio 
of 50, though it is clear that values exceeding this limit are possible 
for different combinations of pulse flux densities and telescopes. 



single pulse searches, using Equation ([5]) as a guide. As 
discussed in ij2.1l detecting a periodic pulse train at the 
distance of M33 is likely to be possible only for pulsars at 
the bright end of the pulsar pseudo-luminosity function. 
In contrast, with their much greater luminosities, giant 
pulses should be detectable to much larger distances than 
M33 with higher sensitivity systems (or at the distance of 
M33 by less sensitive systems). 

Tabic [3] and Figure [5] summarize the detection possibil- 
ities with future instruments. Given the large number of 
facilities becoming available, we do not attempt to summa- 
rize them all here. Rather, we choose several instruments 
as exemplars of future possibilities. While our choice of 
the specific instruments is somewhat arbitrary, there are 
multiple instances in which various instruments have (or 
are expected to have) similar performance. For example, 
the system sensitivities of the Expanded Very Large Ar- 
ray (EVLA) and the Green Bank Telescope (GBT) are 
the same within 5% {Ssys ~ 10 Jy), and the Giant Metre- 
wave Radio Telescope (GMRT; Ssys ~ 11 Jy) will offer a 
similar sensitivity once its new wide-band instrumentation 
(400 MHz) comes online. A future expansion of the Allen 
Telescope Array (ATA-256) has a system sensitivity that 
is within a factor of 2 of these telescopes. 

In the near future, a substantial increase in sensitivity 
(by a factor '-.^2-3) that will be possible through the new 
instrumentation at Arecibo, such as the Mock spectrom- 
eters and the proposed Arecibo Ultimate Pulsar Process- 
ing Instrument, offers significant promises for conducting 
renewed searches. Yet another promise is the proposed 
phased-array feed for Arecibo that will provide 42 beams 
instead of 7 ALFA beams; this will enable a much larger 
instantaneous sky coverage and therefore more efficient 
searches of extended objects such as M33. 

Furthermore, with an enormous increase in affordable 
computing power, it will become feasible to do searches 
with coherent dedispersion. This is especially important 
when searching for giant pulses from Crab-like pulsars, 
which are typically 1-10 /zs in duration at 1.4 GHz, with 
the brighter ones being the na rrower pulses (|Bhat et al.l 
120081: iPopov fc StappersI [20071 ). In such cases, the de- 
tection sensitivity can be substantially boosted through 
matched filtering for very narrow pulses. This presents ex- 
citing, albeit computationally challenging, opportunities if 
there are Crab-like pulsars residing in external galaxies. 

As can be seen from Figure [51 the distance to which 
giant pulses can be detected can be approximately nor- 
malized by Arecibo. Telescopes that have roughly 10% of 
the sensitivity of Arecibo (EVLA, MeerKAT) can probe 
to nearby Galactic satellites, but probably cannot reach 
as far as M33 (or M31); telescopes comparable to Arecibo 
(FAST, SKA Phase 1) can probe essentially the entire Lo- 
cal Group; and a telescope with 10 x the sensitivity of 
Arecibo (SKA Phase 2) can reach well beyond the Local 
Group. 

We have used the nearby galaxy catalog of Karachentsev 
et al. (2004) to assess other potential targets; Figure [TOj 
presents a visual summary. We define a potential target 
to be a galaxy at a distance such that a giant pulse that 
is as bright as 20 kJy at the Crab's distance (~2 kpc) will 
produce a signal-to- noise ratio of at least 5 (Equation [2|). 
In the northern (southern) hemisphere, there are 21 (14) 
galaxies that are potential targets of an Arecibo-class tele- 
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Fig. 10. — Galaxies within 3 Mpc accessible to near-future radio telescopes for searches for giant pulses, based on the nearby galaxy catalog 
of Karachentsev et al. (2004), and assuming a Crab giant pulse of 20 kjy in peak flux density, which will be detected with a S/N of 5 
at M33s distance with Arecibos sensitivity (Equation (2)). In both plots, the azimuth is right ascension and galaxies are projected to the 
celestial equator along lines of constant declination. In both panels, the dashed orange circle shows the distance to which a telescope having 
approximately 10% of the sensitivity of Arecibo (e.g., EVLA, GBT, ATA-256, MeerKAT) could probe and the dot-dash blue circle shows the 
distance to which an Arecibo-class telescope could probe. The full SKA could probe to an even greater distance than shown. {Left) Northern 
Hemisphere; (Right) Southern Hemisphere. 



scope; in both hemispheres, there are 10 galaxies accessible 
to telescopes having approximately 10% of the sensitivity 
of Arecibo (e.g., EVLA, GBT, ATA-256, MeerKAT). Of 
course, not all galaxies in Figure [10] will necessarily host 
giant-pulse emitting pulsars. Nonetheless, there are many 
more targets to consider than just M33. 

Despite the numerous challenges it poses, finding pulsars 
in other galaxies will mark a major breakthrough, given 
the multi-pronged scientific benefits that will ensue. For 
instance, pulsars likely to be detected in such searches will 
be young pulsars with high luminosities, hence their as- 
sociations with any cataloged supernova remnants in the 
host galaxies may yield valuable information about the 
star formation rates as well as the pulsar birth rates and 
populations in those galaxies. Moreover, such extragalac- 
tic pulsars, through their dispersion, scattering and ro- 
tation measures, will also provide information about the 
baryon density and magnetic field in the intervening inter- 
galactic medium. With the sensitivity of the full SKA, it 
may even be possible build a partial census of pulsars in 
many nearby galaxies. 

8. SUMMARY AND CONCLUSIONS 

We have conducted an extensive and systematic search 
for pulsars and transient sources in M33 using Arecibo's 
seven-beam receiver system, ALFA. The galaxy, which is 
located at ~0.95 Mpc and extends over roughly 73' x 45' 
in the sky, was covered with 12 ALFA pointings, or 84 pix- 
els in total, each of which was searched for 2-3 hr dura- 
tions. We searched for both periodic and aperiodic signals, 
spanning up to 1000 pc cm""^ in dispersion measure, and 
over timescales from ~50 /its to several seconds. While no 
signals of periodic or transient nature have emerged from 
these searches, they illustrate several underlying challenges 



and difficulties in carrying out such searches as well as im- 
proved efhciencies achievable through the availability of 
simultaneous multiple beams. 

For periodicity searches, our integration times of 2 hr 
translate to upper bounds of 3-10 fi.Jy in flux densities at 
M33's distance, which means very few pulsars potentially 
detectable even in an optimistic scenario of M33 harboring 
a pulsar population akin to that in the Milky Way. For 
short-duration transients, the limiting peak flux densities 

— 1/2 

from our analysis are S'pk- 20 mJy Wms and the implied 
event rate limits correspond to an upper bound of about 
one Crab-pulsar-like object per kpc'^ in M 33. The limits 
on bright, one-off millisecond bursts (e.g., iLorimer et al.l 
1200711 are not as constraining as those from the PALFA 
survey, however still signifi cant and reinforce the conclu- 
sion of iDeneva et al.l (|2009l) . i.e., the source population is 
unlikely to be homogeneous and isotropic in the sky. 

Despite the numerous technical challenges inherent to 
such searches, the future looks promising. Several new 
radio facilities are currently under construction or in com- 
missioning(e.g., LOFAR, FAST, MeerKAT); their higher 
raw sensitivities (i.e., lower Sgys ) and larger bandwidths 
mean more sensitive searches will be possible. Even the 
existing instruments such as EVLA, GBT and GMRT will 
offer opportunities in the near-term future with the avail- 
ability of new instrumentation which will provide much 
larger bandwidths than currently possible (400-800 MHz). 
Arecibo itself will soon offer renewed opportunities, with a 
substantial increase in sensitivity (by a factor ~2-3) that 
will be made possible through the new Mock spectrome- 
ters and other wideband instruments. 

In the longer term future, instruments such as the Five- 
hundred meter Aperture Spherical Telescope (FAST) and 
SKA Phase 1 will bring exciting opportunities for even 
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more sensitive searches. While their sensitivities may be 
comparable to that of Arecibo, their sky coverages will be 
much larger and as a result they will essentially probe the 
entire Local group. Furthermore, with impressive techno- 
logical advances being made in affordable supercomputing, 
searching with coherent dedispersion will eventually be- 
come feasible, which will substantially boost the detection 
sensitivity for giant pulses from Crab-like pulsars that may 
reside in the galaxies. Discoveries of extragalactic pulsars 
and transients will undoubtedly yield a multitude of sci- 
entific benefits, including understanding the populations 
of such objects in the host galaxies as well as serving as 
useful probes of the intervening intergalactic medium. 
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Note. — The coordinates given for each pointing are for the central ALFA beam. 



Table 2 

Statistics of RFI and candidate events 



Quantity 


Matched filtering Time domain clustering 
Short events Long events 




RFI events (classification 4) 


No. of data segments 

No. of events 

% of observing time 


143 40 153 
207 83 375 
0.06 0.024 0.11 




Interesting events (classification 1 or 2) 


No. of data segments 

No. of events 

% of observing time 


72 16 23 
124 18 27 
0.035 0.005 0.008 



Note. The classification into "short" and "long" is based on the event timcscale 
being ;S 10 ms and k. 10 ms, respectively. "Data segment" corresponds to a 4-s chunk 
around the event time; often multiple events are detected within a single 4-s chunk, 
hence the total number of events larger than the number of data segments. 
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Table 3 

Future Searches for Extragalactic Pulses 



Telescope 


V 












(Jy) 
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(at Dm33 ) 
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3.2 
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11.0 


400 


2.7 


FAST 
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400 


13.6 


MeerKAT 
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1.9 


SKA Phase 1 


1.4 GHz 


3.6 


500 


9.3 


SKA Phase 2 


1.4 GHz 


0.3 


1000 
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Note. — The 


instruments listed 


are the 


Arecibo Ob- 



servatory, the Low Frequency Array (LOFAR), the Ex- 
panded Very Large Array (EVLA), the Giant Metrewave 
Radio Telescope (GMRT), the Five- hundred- meter Aper- 
ture Spherical Telescope (FAST), the Karoo Array Tele- 
scope (MeerKAT), the SKA Phase 1 (SKA-1), and the SKA 
Phase 2 (SKA-2). As many of the telescopes are still under 
construction or have just started commissioning, parame- 
ters listed should be considered indicative but not defini- 
tive. {S/N)gp denotes the predicted signal-to-noise ratio 
(at M33s distance, DM33) for a giant pulse that is as bright 
as ~20 kJy at the Crabs distance (~2 kpc). 



